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We consider quasi one dimensional spin-1 Heisenberg chains with crystal field anisotropy in a 
uniform magnetic field. We determine the dynamical structure factor in various limits and obtain a 
fairly complete qualitative picture of how it changes with the applied field. In particular, we discuss 
how the width of the higher energy single magnon modes depends on the field. We consider the 
effects of a weak interchain coupling. We discuss the relevance of our results for recent neutron 
scattering experiments on the quasi-lD Haldane-gap compound NDMAP. 

I. INTRODUCTION 

Recent years have seen a resurgence of interest in field-induced "magnon condensation" in gapped quasi one- 
dimensional quantum antifcrromagnets. In particular a series of ESR and neutron scattering experiments have been 
carried out on the Haldane-gap 1 chain compounds NDMAP 2-7 and NDMAZ 8 . The main motivation for these ex- 
perimental studies is the observation that a spin-1 Heisenberg chain undergoes a quantum phase transition between 
a gapped spin-liquid phase and a gapless Luttinger liquid phase at some critical value H c of the applied magnetic 
field H 9 ' 10 . The ground state of the spin-1 Heisenberg chain is a spin singlet and excitations are described in terms 
of a gapped S = 1 triplet of magnons. When a magnetic field is applied, the triplet splits due to the Zeeman effect 
and one of the magnon gaps is driven to zero at H c . For H > H c the ground state is magnetized and excitations 
are gapless. If interactions between the magnons were absent, the transition at H c could be understood as a Bose- 
Einstein condensation of magnons. In the spin-1 Heisenberg chain there is an interaction between magnons, which 
fundamentally changes the ground state for H > H c from a condensate of bosonic magnons to a Luttinger liquid, 
which can be regarded as a one dimensional version of an interacting Bose condensate. The transition at H = H c 
is in the universality class of the commensurate-incommensurate (C-IC) phase transition 11 . The magnetic response 
of the isotropic spin-1 chain in strong fields H > H c has been analyzed in some detail in Refs [12-14]. In Appendix 
A we use the nonlinear sigma model description of the isotropic spin-1 chain to derive explicit expressions for the 
dynamical response functions in the low- field phase H < H c . In many S = 1 compounds such as NENP, NDMAP 
and NDMAZ strong crystal field anisotropies are present, which lead to a zero-field splitting of the magnon triplet 
comparable in magnitude to the Haldane gap itself. These anisotropy effects lead to more complex behavior and a 
richer phase diagram. The purpose of this work is to determine the dynamical response of Haldane-gap compounds 
in the presence of such crystal field anisotropies. The relevant lattice Hamiltonian is of the form 

H = j'YlSj- S j+1 - H Sj 

+ E[(S?) 2 - (SV) 2 ], (1.1) 



where < — E < D. In zero field there are three magnon modes in the vicinity of the antiferromagnetic wave number 
q = jL. with gaps Ai < A 2 < A3. For simplicity we will mainly concentrate on the case where the magnetic field is 
applied along the z-axis, i.e. 

H = He z . (1.2) 

The lattice Hamiltonian (1.1) (with field applied along the z-direction) exhibits two discrete symmetries which play 
an important role in the following: 

1. Rotation by ir around the z-axis (R4): 

sf -> -sf , sy -> -sy , s z -> s*. (1.3) 



2. Translation by one site (Tr): 



Sf — Sf +1 , a = x,y,z . (1.4) 
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The model (1.1) is difficult to analyze directly by analytical methods. However, progress can be made by concen- 
trating on the low-energy regime, which can be studied by means of semi-phenomcnological descriptions in terms of 
continuum models. Two such models have been used in particular, namely (i) a bosonic Landau-Ginzburg model 10,14 
and (ii) a theory of three coupled Majorana fcrmions 15 . In the present work we go beyond the original works 10,14 ' 15 
by (1) discussing the decay of high energy magnon modes, (2) applying methods of intcgrable quantum field theory 
to the discussion of structure functions and (3) taking into account the effects of inter-chain interactions. 

The outline of this paper is as follows. We start by reviewing known results on the spectrum of the Majorana 
fermion model in section II B. We then investigate the role played by interactions in sections II C and II D and in 
particular show that interactions generate a finite lifetime for one of the magnon modes. In section III we derive 
analogous results in the framework of the Landau-Ginzburg theory. We then turn to a more detailed analysis of the 
magnetic response functions in the low-energy regime in the vicinity of the quantum critical point at H c in section IV. 
Section V gives a detailed account of the low-energy regime in the high-field phase for small crystal-field anisotropy. 
The effects of interchain coupling are investigated in section VI and a summary and discussion of our various results 
is given in section VII. A variety of technical details are presented in Appendices A-E. 



II. MAJORANA FERMION MODEL 



In Ref.[15] A.M. Tsvelik proposed a description of the spin-1 Heisenberg chain in terms of a field theory of three 
right and left moving Majorana (real) fermions R a = R\, L a = L\, a =1,2, 3. The Hamiltonian is given by 
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H = - ^2 v [Lad x L a — R a d x R a ] — A a [R a L a — L a R a ] 



o=l 

+ i_ J- s abc H a (L b L c + R b R c ) + H' , (2.1) 

where H' describes anisotropic current-current interactions 

H' = Y,9aJ a J a ■ (2-2) 

a 

Here the currents are bilinears in the Majorana fermions 

J a {x) = - l -e abc [L b L c + RbR c ] ■ (2.3) 

We will take the magnetic field along the 3-axis and concentrate on the case 

Ai < A 2 < A 3 . (2.4) 

The generalization of our results to other situations is straightforward. The smooth and staggered components of the 
spin operators are defined by the decomposition 

Sf^ J a (x) + (-iyn a (x) , (2.5) 

where x — jao- The smooth components are equal to the currents, where a = 1,2,3 correspond to a = x, y, z. The 
staggered components are expressed in terms of Ising order and disorder operators 

n x (x) oc a 1 (x)fi 2 (x)fi 3 (x) , 
n v (x) oc ii 1 (x)a 2 (x)ii 3 (x) , 

n z (x) oc fi 1 (x)fj 2 (x)a 3 (x) . (2.6) 
We note that the signs of the mass terms in (2.1) are such that in zero field 

(a a (x)}=0, (v a (x))^0. (2.7) 
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A. Symmetries 

The Hamiltonian (2.1) inherits the discrete symmetries Tb and R* from the lattice model. The latter is realized 



as 



R4 : R a 



-R„ 



— (7 a , fi a — > n a , a = 1,2. 
The translation symmetry by one site turns into a discrete Z 2 symmetry in the continuum limit 

J a {x) ► J a {x) , n a {x) ► -n a {x) , 

and may be realized as 

Ra y Ra 7 L a > L a , 

a a — > -a a , fi a — > n a , a = 1, 2, 3. 
It is convenient to combine this symmetry with R4 into the following Z 2 symmetry 



T fl R* : R 3 

3 



a — > — (T 

The full symmetry of the Majorana model (2.1) is thus Z 2 ® Z 2 . 



— -R3 ) L 3 — > — L 3 
.3 



M 3 > A* 3 • 



(2.8) 
(2.9) 

(2.10) 

(2.11) 



B. Spectrum in the absence of interactions 

In the absence of interactions, i.e. g a — in (2.2) the Hamiltonian can be diagonalized by means of a Bogoliubov 
transformation 15 . In the following we review some relevant formulas. As the magnetic field is along the 3-direction 
only the first and second Majoranas couple to the magnetic field. The third Majorana gives rise to a fermionic 
single-particle mode with dispersion 



cc 3 (fc) = ^/A 2 3 +v 2 k 2 . 
The first and second Majoranas are conveniently combined into a complex fermion \& 



Ri = 
Li = 



V2 



R2 = 



iV2 



Lo = 



V2 ' iV2 ' 

The Hamiltonian density describing the first and second Majoranas takes the form 



n 



12 



— IV 

H r 

—im 



iA 



#j,# L - /i.e. 



- h.c. 



where 



Introducing a mode expansion 



. A 2 - Ai A 2 + Ai 
A = , m = . 



*«(a;) 
V L (x) 



= J o ^-[e ik *a(k)+e- ik */3Hk)] , 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



(2.16) 
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we may express the Hamiltonian density (2.14) as 
where 



r°° dk* t 

Jo ^ ^ 7a 

a, 6=1 



(fc)M ab7b (fc) 



la (k) = (a(k),a^(-k),p(k),^(-k)) a , 
M = 



vk + H 


iA 





—im 


-iA 


-vk-H 


im 








—im 


vk-H 


iA 


im 





-iA 


-vk + H 



Now we perform a Bogoliubov transformation with a unitary matrix U(k) 

U(k) 



( c+{k) \ 
cl(-k) 
c-(k) 
\cl(-k)J 



at(-jfe) 
(3(k) 
\PH-k). 



to bring H\i to a diagonal (normal ordered) form 

^5>a(*0 cl(k) Ca (k), 



-°° — a =± 

m 2 + A 2 + H 2 +v 2 k 2 



±2 v / m 2 A 2 + H 2 (m 2 + v 2 k 2 ) . 
The gap w+(0) increases monotonically with H, whereas w_(0) decreases and vanishes at a critical field 

H c = ^A~ 2 - Vm 2 - A 2 . 
The corresponding critical point is in the Ising universality class. For fields H < H C 2 



H c2 = m 



I— 1 


I— 1 


A 2 


2^ 




2 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



the minimum of the dispersion occurs at k = 0. In the vicinity of H c and at small momenta the dispersion is 

approximated as 



Z(fc) « ^(H-H c ) 2 



m' 



+ v' 



A 2 H c (l 



m^ 



+ 



m"- 



■(H c -H) 



k 2 . 



Hence the gap vanishes linearly with H — H c in agreement with Ising critical behaviour 



A 2 

_(0)^(H-H c )Jl-—. 

V m z 

In order to see the Ising criticality described by (2.23) the magnetic field must be sufficiently close to H c 

A 2 



H — H r < 



(2.23) 



(2.24) 



(2.25) 
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As expected this scale is set by the anisotropy A. For H > H c2 there are two degenerate minima of w_(fc) at some 
incommensurate wave numbers ±kp with 



For k w ±kp we have 



vki 



H 2 



m 2 A 2 



m 



tf 2 



uj-(k F ) = A\ 1 - 



if 2 ' 



L0 2 _(k) 



1 



TYL 
IP 



+ 



v v 2 ki 



-i 2 



H 



k F f. 



(2.26) 



(2.27) 



As was recently pointed out by Wang 26 , these results for the dispersions suggest that a cross-over between Ising and 
C-IC critical behaviour occurs as a function of H . For H very close to H c wc encounter Ising critical behaviour, which 
crosses over to C-IC behaviour for H > H c2 . 



C. Interactions: self-consistent mean-field treatment 



So far we have neglected the four-fermion interactions (2.2) altogether. As a first step of taking interactions into 
account we may treat them in a self-consistent mean- field approximation (SCMF). The following expectation values 
are compatible with the discrete Z 2 symmetries 

(L a R a ) ? , (LiL 2 ) ^ , (i?!i? 2 ) ± , 

(JiiL 2 ) ? , (R 2 L 1 ) + . (2.28) 

Decoupling the four fermion terms leads to a renormalization of the mass parameters 

Ai — » Ai = Ax- 2ig 2 (L 3 R 3 ) -2ig 3 (L 2 R 2 }, 
A 2 — » A 2 = A 2 -2i.g 1 ( J L 3 i?3>-2«.g 3 (Lii 1 , i), 

A 3 — » A 3 - A 3 - 2^.g 3 ( J L 2 i? 2 ) - 2ig 2 (L 1 R 1 ) . (2.29) 
The magnetic field terms are changed to 

# L Li£ 2 + iffl i?ii? 2 , (2.30) 

where = iJ + 2ig 3 (RiR 2 ). and i?_R = H + 2ig 3 {L\L 2 ) . Finally, two new terms are generated 

»Ai RiL 2 + i\ 2 LiRa , (2.31) 

where iX\ = 2g 3 {L\R 2 ) and i\ 2 = 2g 3 (R\L 2 ). The resulting mean-field Hamiltonian is quadratic in the Fermi fields 
and can again be diagonalized by a Bogoliubov transformation. Let us denote the ground state energy obtained in 
this way by Eqs. The expectation values are determined self-consistently, e.g. 

—i(R a L a ) = ^ . (2.32) 

dA a 

The SCMF procedure is easily implemented once the appropriate couplings g a are specified. However, in order to 
keep matters as simple as possible we will assume from now on that the g a are small and as a result the differences 
between the free theory and the SCMF theory are negligible. The main qualitative effect of nonzero g a is to make the 
gap of the third Majorana magnetic field dependent. Such a dependence is observed for example in experiments on 
NDMAP 4 ' 5 and implies the presence of interactions in the framework of the Majorana fermion model. 



D. Decay Processes in the Low-Field Phase 



Within the SCMF approximation the role of the current-current interactions is merely to induce slight changes 
of the dispersion relations of the three coherent single-particle magnon modes. As we will now show, treating the 
interactions beyond the SCMF approximation leads to the damping of one of the magnons. 
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The analysis of the spectrum summarized above establishes that there are three different types of magnons, which 
we will refer to as M 3 , M + and M_ respectively. The corresponding dispersion relations are w 3 (k) (2.12) and 0J±(k) 
(2.20). The interaction of these modes is described by the term Tt' in the Hamiltonian (2.1) and involves four particles. 
As u)-(k) can become very small as the magnetic field is increased from zero, the decays M 3 — ► M_Af_Af_ and 
M + — ► M_M_M_ become kinematically allowed for sufficiently large magnetic fields. However, the decay of M 3 is 
forbidden by the symmetry T#R*: M 3 is odd under this symmetry whereas M± are even. Essentially we are using 
the fact that all 3 magnon modes only exist near wave-vector so that they can only decay into an odd number of 
magnons. The decay of M 3 into any odd number of M_'s would be inconsistent with the R% spin rotation symmetry. 
On the other hand the decay M + — > M_M_M_ is allowed by symmetry. The process becomes kinematically possible 
as soon as the magnetic field exceeds a critical value which is defined by 

w+(0) = 3w_(0) . (2.33) 

Solving for we find 



H d = \j— -A2 . (2.34) 

As long as 2A < m the decay process will occur for H > and from now on we will assume that this is the case. 

We note that in zero field there are no decay processes even if the gaps are such that they are kinematically allowed 
(3Ai < A 2 ). The reason is that for H = the lattice Hamiltonian (1.1) has additional spin rotational symmetries 
around the x and y axes by 180 degrees. In combination with the translation symmetry these induce symmetries of 
the form (2.11) for the Majoranas 1 and 2 individually rather than the combination (2.8). This additional symmetry 
forbids decay processes. 

Inserting the mode expansions (2.16) into the expression for H' , integrating over the spatial coordinate and then 
carrying out the Bogoliubov transformation (2.19) generates several terms quartic in the fcrmionic creation annihila- 
tion operators c±(k), c±(k). The most interesting one describes the decay of a M+ mode into three M_ modes and 
is of the form 



t/ f°° dk idk 2 dk 3 

V = 93 (2n)* A* 1 ' *»■*») 



(2tt) 3 

xcl{k 1 )cl(k 2 )cl(k 3 )c + (k 1 +k 2 + k 3 ). (2.35) 

Here / is an antisymmetric function of its arguments and at small momenta is of the form 

f(k u hi, k 3 ) ~ C (fci - hi){k x - k 3 )(k 2 - k 3 ). (2.36) 

The constant C is a complicated function of A, m and H . The differential rate for the decay of a M + particle with 
momentum p into three M_ particles with momenta Pi,p 3 ,p 3 is 

3 3 

dT = 2 7 r|M| 2 5(p-^p j )^+(p)-E w -fe)) 
j=i i=i 

dpidp 2 dp 3 



3! 



(2.37) 



where the factor of 3! is introduced to account for the fact that the three particles in the final state are indistinguishable. 
In the Born approximation the transition matrix element M is 



3 3 

M = ^(o\f[c-(p 3 )v ciipmfa-^Tp^y 1 

f(pi,P2,p 3 ). (2.38) 



2tt 

Taking the M + particle to be at rest, i.e. setting p = 0, we obtain 

r.2 



60 f 

r = / d Pl d P2\f(Pl,P2,-Pl - P2)\ 2 

x <y(w + (0)-w_(pi)-w_(p2)-o;-(pi+P2)). (2.39) 
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In order to simplify matters further, we concentrate on the case where the magnetic field is close the critical field H d 
at which the decay M + — ► M_M_M_ first becomes kinematically possible. In this regime we have 



w+(0) - 3u_(0) < w_(0) . 



(2.40) 



Then the momenta in (2.39) have to be small in order to satisfy the delta-function and we may use the expansions 
(2.36) for / and 



» w w_(0)+ap 2 +O(p 4 ) 
H 2 



a = v 



1 - 



[^-(O)]" 1 . 



mVH 2 + A 2 _ 

This leads to the following expression for the decay rate in the regime H > H d , ^-1<1 



(2.41) 



6gj C 2 4tt 
2^ (25) 4 ^3 

ff 2 V3C : 



[w + (0)-3w_(0)] 3 



"2" 


4 




A 2 " 


a 




1 - 


Am 2 



(2.42) 



We find that the decay rate is proportional to (H — H d ) 3 and is therefore quite small in the vicinity of H d . 



III. LANDAU-GINZBURG (LG) MODEL 

A different approach to studying the spin-1 Heisenberg chain with crystal field anisotropics in a magnetic field 
was used in Refs [10,14]. It is based on the nonlinear sigma model description of the spin-S Heisenberg chain in 
terms of the three-component field (p describing the staggered components of the spin operators and the subsequent 
approximation of replacing the constraint Cp 2 = 1 by a |</?| 4 interaction. The LG Lagrangian density is 10,14 

2v \ dt W 2 V dx 



-E 



A yl-\\0\ 4 - (3.i) 



a=l 

We again take the magnetic field to point along the 3-axis 

H = He 3 . (3.2) 

It then follows from (3.1) that <p3 couples to the magnetic field only via the \\<f\ 4 interaction. The Landau-Ginzburg 
theory inherits the discrete symmetries (1.3) and (1.4) from the underlying lattice model. As is the Majorana fermion 
model it is convenient to combine the translational symmetry by one site T^j with the rotation around the z-axis by 
7r R4 and obtain the following two Z 2 symmetries 

R*: Va^-<Pa, a -1,2. (3.3) 

T fl R4 : p 3 — » -<p 3 . (3.4) 
The resulting Z 2 (8> Z 2 symmetry is the same as for the Majorana fermion model. 

A. Spectrum and Mode Expansion 

Neglecting the \<p\ interaction the spectrum can be determined by solving the classical equations of motion 



>: ' -\ a+ 2**°H b °£=0. (3.5) 



W 2 - v2 ^ +a °- h2 



dt 
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One finds that there are three magnon modes M 3 and M± with the following dispersion relations 10,14 



UJ 3 (k) = X /Al+«2fc2 



oj±(k) 



A 2 4- A 2 
Ht+ ^Z 2 +v 2 k 2 



±\ 2H 2 (Al + A 2 + 2v 2 k 2 ) + 



A? - A 2X 1 



(3.6) 



The low energy M_ mode with dispersion u)-(k) becomes gapless at a critical field H c = Ai. The resulting critical 
point is in the universality class of the two-dimensional Ising model 14 . 
The scalar fields ipi^ have the following mode expansions 



ipi(t,x) = Yl I dk 

a=± J 

<P2(t,x) = ^2 / dk 

a=± J 



A la (k)e- lu ^ t+lkx a a (k) + h.c. 

ii 2n (^- M °W' +ih a tt (fc) + h.c. 
y/4nu a (k)/v 



where A* aa (k) = A aa {k). Here a and a) obey canonical commutation relations 

[a a (k),ai(p)} = 5 at 3 S(k - p) , 



(3.7) 



(3.8) 



and the amplitudes A a ±(k) are fixed by the requirements that (i) the fields (p a fulfil the equations of motion (3.5) and 
(ii) the fields ip a and the conjugate momenta Il a = \ ^§f + (H x tp) a fulfil canonical commutation relations 



[<p a (t, x),<p b (t, y)} = 0, [n a (t, x), n b (t, y)]=Q, 
[fa(t, x),U b (t, y)] = iS ab S(x - y) . 



(3.9) 



We find 



A\ + (k) 



3H 2 + A{ +v 2 k 2 - w_(fc) 2 
oj+(k) 2 -w_(fc) 2 



AUk) = l-Al +{ k) = ^?-^-*l-v 2 k 2 



A 2 g(k) 

A la {k) 



u+(k) 2 -u_(k) 2 
cu a (k) 2 + H 2 - A 2 - v 2 k 2 



2Hcj a (k) 



, a = ±. 



(3.10) 



Eqns (3.10) allow us to deduce the polarizations of the modes corresponding to the dispersion ui±(k). For example, 
at H = the M_ mode is polarized entirely along the 1 direction and the M + mode along the 2 direction. On the 
other hand, as H — > Ai wc find that 



A 2 -(0) 
A 2 _(0) 



0, A 2 + (0)^1, 



A 2 



3A 2 + A 



5, ^+(0) 



4A 2 



3Aj + Aj. 



2 • 



(3.11) 



Hence tp 2 couples only to the M + mode whereas ip\ couples to the M + mode as well as to the M_ mode with a 
strength set by the anisotropy. 



B. Decay Processes 

In the absence of the nonlinear |<p| 4 -term the LG model describes three coherent magnons M 3 , M± with corre- 
sponding dispersion relations (3.6). Inclusion of the \0\ 4 term generates interaction terms involving four particles. As 
was the case in the Majorana fermion model, ui—(k) can become very small when the magnetic field is increased and 
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as a result the decays M3 — ► M_M_M_ and M + — > M_M_M_ become kinematically allowed. The decay of M3 
is forbidden by the symmetry T^R^ (3.4): M 3 is odd under this symmetry whereas M± are even. On the other hand 
the decay process M + — » 3Af_ is allowed when the magnetic field is larger than 



-[ 
16 L 



[A? + A 2 ] ± A^i 3[A 4 + a4] + ioa?a; 



The interaction describing this decay process is given by 

Aw 2 r 00 

V = — — / dk 1 dk 2 dk 3 g(k l7 k 2 ,k 3 ) 

xal(ki)aL(k2)aL(k 3 )a + (ki + k 2 + fc 3 ), 
where g(k\, k 2 , k%) is a symmetric function of its arguments. Its zero momentum limit is 



(3.12) 



(3.13) 



5(0,0,0) 



A!_(0) 3 A 1+ (0) 



UJ 3 _UJ + 



lo 2 _+H 2 -A 2 ^ 2 
2Hlu_ 



u 2 _+H 2 -A( u 2 +H 2 -A{ 

1 H — - 

2Huj- 2Hlo + 



C 



UJ 3 _UJ + 



= .(3.14) 



where lo± == u ± (0). As a first approximation we neglect all interactions except V and calculate the decay rate 
M + — > 3M_ in the Born approximation. 

The differential rate for the decay of a M + magnon with momentum p into three M_ magnons with momenta 
P11P31P3 is again given by (2.37), where 



M = (0\Ha-( Pj )Val(p)\0) Up-^p,-) 

3=1 V •J= i 

3Xv 2 , 
= g{PuP2,P3)- 



(3.15) 



Taking the M + magnon to be at rest, i.e. setting p = 0, we obtain 

3X 2 v 4 f 

T = — — / dp 1 dp 2 \g(pi,p 2 ,-pi -P2)\ 2 

x S(u+(0) - u-(pi) - w_(p 2 ) - +P2))- 
For H slightly larger than H d we again have 

w+(0) -3w_(0) < w_(0) , 



(3.16) 



(3.17) 



and the momenta pi, 2 in (3.16) have to be small in order to satisfy the delta-function. We may use the expansions 
(3.14) for g and 



w_(p) wco_(0)+V+O(p 4 ) , 



a = 



2w_(0) 



2# 2 



2 J ff 2 (Af + A2) + (^i^)2 



(3.18) 



This leads to the following expression for the decay rate in the regime H > H d , jf — 1<1 

~ cm>_(0) 3 w + (0)' 



(3.19) 



The result (3.19) would suggest that the decay rate switches on suddenly at a finite value as soon as H becomes 
larger than H d . The underlying reason for this jump is the restricted phase space of a one-dimensional system. The 
"free boson" result (3.19) for the decay rate is dramatically different from the result obtained in the framework of 
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the Majorana fermion theory. This poses the question whether (3.19) is robust if we take into account interactions 
among the magnons in the final state. This amounts to resumming the leading infrared divergences in a perturbative 
expansion in A. If we assume that in the low-energy limit the M_ degrees of freedom in the Landau-Ginzburg 
model can still be mapped onto a Bose gas with (5-function interactions 16 , the result of a such a resummation can be 
determined by exploiting the fact the wave functions of the M_ modes reduce to a free fermion form in the limit of 
small momenta 19 . A three-particle state in the position representation can be written as 



|Ai,A 2 ,A 3 ) = J dx 1 dx 2 dx 3 x(xi,x 2 ,x 3 ) 

x al(xi)al(x 2 )al{x 3 )\0), (3.20) 



where the wave-function is given by 



X3(xi,x 2 ,x 3 ) = 3 /23 , Yl sgn( Xj -x k ) 
' ' j<k=l 



sgn(P)e^=i <Ap ^'. (3.21) 



x 
Pes 3 



Here P denotes a permutation of three elements and 5*3 the symmetric group of degree 3. In momentum space we 
have 



3 

|Ai, A2, A3) - / dfc iy 3 fl J^±- a t_ (Al _ k2 _ fc 3 )a t _(A 2 -k 1+ fc 3 )a t _(A 3 + h + fc 2 )|0). (3.22) 
Taking (3.22) as the final state, the matrix element of the decay vertex is 

= 6{p _ g 3 =i Aj) = — J H ^JT7 2 g(Ai - ft - A 2 - gl + ft, A 3 + gl + ft ) . (3.23) 

It follows from the definition (3.23) that the matrix element is an antisymmetric function of Ai, A 2 , A 3 . As long as we 
are interested in the decay rate for H close to i?<a, we may expand M for small Xj. The leading term antisymmetric 
in Xi 2.3 is then 

C"(Ai - A 2 )(Ai - A 3 )(A 2 - A 3 ) , (3.24) 
For small Xj the matrix element is thus equal to 

M = ^!rr (Aj .-A fe ). (3.25) 

j<k 

Following through the same steps as before, the decay rate for H close to with the M + magnon initially at rest is 
found to be 



V3X 2 C' 2 



4 



4a 4 



[u)+{0) -3w_(0)f cx [H-H d y\ (3.26) 



Comparing the decay rate (3.26) to the result (3.19) we see that interactions in the final state have a dramatic effect: 
rather than turning on at a finite value as soon as H exceeds H^, the decay rate (3.26) actually vanishes at and 
exhibits the same power law behavior for H > as the decay rate in the Majorana fermion model. This may suggest 
that the dependence of the decay rate on H — H d is a robust result that holds for the underlying lattice model as well. 



IV. VICINITY OF THE ISING CRITICAL POINT 



As we have seen above, both the Majorana fermion model and the Landau-Ginzburg theory lead to an Ising critical 
point at some values H c of the applied magnetic field. In the simplest approximations where interactions are neglected, 
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4 




FIG. 1: Magnon gaps as functions of the applied field for the Majorana fermion model with g a — 0. The gaps in zero field 
are chosen as Ai = 0.8m, A2 = 1.2m and A3 = 3.5m. The square indicates the window in energies and fields in which the 
description in terms of an effective Ising model is appropriate. 



the actual values for ff c are rather different, but as we are interested only in robust features this numerical difference 
does not really concern us here. More importantly both theories predict that the low-energy behaviour in the vicinity 
of the critical point is described by an off-critical Ising model with Hamiltonian 

xv 

H = — [Ld x L - Rd x R] - im RL . (4.1) 

Here the mass parameter that parametrizes the deviation from criticality is equal to the smallest magnon gap 

\m \ ~ uo-{k = 0,if) . (4.2) 

The description by (4.1) is appropriate at low energies u> <C min{cj + (fc = 0, if), u>3(k — 0, if)} and if sufficiently 
smaller than the critical field if C 2 (2.22) at which incommensurabilities develop. The Hamiltonian (4.1) exhibits a Z2 
symmetry 

R — >-R, L — >-L, (4.3) 
under which the Ising order (<r) and disorder (/1) parameters transform as 

H <H C : a — ► -a , 

M — ► M i 
H > H c : a — > a , 

M — -M ■ (4.4) 

In order to determine the magnetic low-energy response (in the vicinity of the antiferromagnetic wave number) we 
need to express the staggered magnetizations n a in terms of operators related to the Ising model (4.1). As was pointed 
out in Ref.[14], the dominant contribution to the staggered magnetization in the x-direction should be the Ising order 
parameter field 

n x =C X (H) a + ... (4.5) 

Here C X (H) is an unknown constant and the dots indicate contributions from less relevant operators. The underlying 
reason for the identification (4.5) is simply that in the ordered phase of the Ising model (which corresponds to if > if c ) 
one has a nonzero expectation value (a) ^ 0, whereas in the disordered phase (which occurs at if < if c ) one has 
(er) = 0. The staggered magnetization in both the Majorana fermion and LG models exhibits the same kind of 
behaviour and it is then natural to expect the identification (4.5). In Appendix D we present arguments that suggest 
that the staggered magnetization along the y-direction has the following low-energy projection 

n y = C y (H) d T a + . . . (4.6) 
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In what follows we chose a short-distance normalization for the field a in the theory (4.1) such that for t 2 + x 2 /v 2 — > 

(a(r, x) a(0, 0)) — ^— L-_. . (4.7) 

We will use the integrability of the theory (4.1) to determine the dynamical structure factor S xx (ui, + q) at low 
energies, where the staggered component of S x is given by (4.5). The yy component then follows from (4.6) 

2 

S yy (u, — +q)oc ^^(u, — +q). (4.8) 



1. H < H c : Low-Field Phase 



By the symmetry (4. 3), (4. 4) only intermediate states with an odd number of magnons contribute to the two- 
point correlation functions of the Ising order parameter field a and hence the staggered structure factor S xx (u), ^ +<?)) 
at low energies. The leading contributions are 20-23 

7T . VA 



S™(cu, - + q)= 6 [w Jm 2 + v 2 q 2 

a o \Jm z + v z q z \ v 

2vA f z ° (tanh(z)tanh(2±^)tanh(2-£)) 2 
' dz 



,2 



" aJo ^/[a; 2 -l-4cosh 2 z] 2 -16cosh 2 z 
+ contributions from 5,7,. . . magnons, (4.9) 

where 

A = C%{H)2h-iA 3 ml , (4.10) 

^2 _ ' 



v 2 q 2 



z = arccosh( — - — ) , 



y 

Here 



/ ^ z — \ — 4 cosh 2 z \ 

arccosh . (4.11) 

V 4 cosh z 



A = 1.28242712910062... (4.12) 



is Glaisher's constant. The three-magnon contribution is always very small. In the frequency interval [0,30too] 
roughly 100 times more spectral weight sits in the single-magnon contribution than in the three-particle one. Hence 
the magnetic response below energies of the order of tens of the magnon gap mo is dominated by the coherent magnon 
contribution. However, if H becomes very close to H c the magnon gaps m a tends to zero. If we are interested in the 
magnetic response at a low (compared to the gap of the second coherent magnon mode) but fixed energy we have to 
take the contributions of intermediate states with 5,7,9,... magnons into account in order to get an accurate result for 

2. H = H c : At Criticality 
At criticality the structure factor exhibits a power-law behaviour 24 

S xx (u, — + q) = Im\B[v 2 q 2 -(Lu + ie) 2 ]~ 7 *}, (4.13) 

where 
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3. H > H c : High-field Phase 

In this phase there is a nonzero staggered magnetization and correspondingly a nonzero expectation value for the 
Ising order parameter 

(a) ? 0, (4.15) 

which results in a Bragg peak for momentum transfer along the chain direction. By virtue of the Z 2 symmetry 
(4.4) only intermediate states with an even number of magnons contribute to S' xx (a;, + q). The leading contribution 
to the inelastic neutron scattering cross section comes from intermediate states involving two magnons 20-22 , which 
leads to the following result 



nsx , 7r vA J s 2 - 4m 2 , 

S™{w, — + q) = 3 -9(s - 2m ) 

+ contributions from 4,6,. . . magnons . (4.16) 
Here A is given by (4.10) and s 2 = J 1 - v 2 q 2 . 

V. THE HIGH-FIELD PHASE FOR WEAK ANISOTROPY 

In general it is difficult to determine the effects of magnon-interactions in quantitative detail. An exception is the 
case of a small anisotropy of the zero-field gaps 

A«m, (5.1) 

where 2m = Ai + A 2 and A = A 2 — Ai. As we will show, if the magnetic field is sufficiently larger than the critical 
field H c and 

A < H - H c < J , (5.2) 

it is possible to determine the interaction effects on the magnetic response at low energies in some detail. More 
precisely, we consider magnetic fields sufficiently larger than the field H c2 (see e.g. (2.22)), at which incommensura- 
bilities begin to develop. The scale defining the low-energy region is the difference A. As is shown in Appendix E, 
the low-energy effective Hamiltonian is given by a sine-Gordon model (SGM) 14 

H = -J- [(^6) 2 + (d x <t>) 2 ] - 2 M cos/36 . (5.3) 

167T 

Here v is the spin velocity, fj, oc A and [3 is a function of the applied magnetic field H . The high-energy cutoff of the 
theory (5.3) is H — H c . At low energies the dominant Fourier component of the transverse spin operators is at q = ^ 

Sf (-l)^cxp(±^e) . (5.4) 

The identification (5.4) is in accordance with the fact that in the high-field phase there is Neel order along the 
x-direction 

((-l)"^)a(cos(^e))^0. (5.5) 
We choose a short-distance normalization such that for \x — y\ — > 

The amplitude A in (5.4) is nonuniversal and not known in general. However, very close to H c 

A < H-H c < H c , (5.7) 

it is given by (see Appendix F) 

A = A'\ a ^(H-H c )V , (5.8) 
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where A' is a field independent numerical constant. The dependence of (5.8) on H — H c is a universal feature of the 
C-IC transition. 

We note that the sign of the cos-term in (5.3) is quite important. Flipping the sign corresponds to a shift O — ► 
+ n/(3, which essentially leads to an exchange of the x and y component of the spin operators in (5.4). 

The value of the parameter j3 is of crucial importance. For the isotropic case (A 2 = Ai) it has recently been 
determined 13 in the framework of the nonlinear sigma model description of the spin-S Heisenberg chain. In the 
isotropic case the high-field phase at H > H c is a Luttinger liquid and (5 is related to the Luttinger liquid parameter. 
It was found that 13 



V2S R (9 F ) ' 



where Sr(0) fulfils the integral equation 

Sr(9) = 1 + 



J-t 



dff S R (9') 



1 



-e F ^ + {9- 0- 

Here 9 F is determined as a function of the magnetic field H by 

r 9 F i 



/\2 



(5.9) 



(5.10) 



d9't{9')— 9 
e(6 F ) = 0. 



(5.11) 



Similarly one may determine the spin velocity 

0.8 



CO- 0.6 




FIG. 2: Parameter (3 as function of the applied magnetic field. 



de(6) 



2np{6) 89 



e=e F 



where p{9) fulfils the integral equation 



p(0) = _cosh(0) + 



(5.12) 



(5.13) 



The parameter (3 as well as the velocity v entering the sine-Gordon Hamiltonian (5.3) may be estimated with a 
good degree of accuracy from their respective values in the isotropic case as long as A2 A ' Al <C 1. On the other hand, 
the agreement between the Luttinger liquid parameter calculated in the nonlinear sigma model and the one of the 
isotropic spin-1 Heisenberg chain in a magnetic field, which is known approximately from DMRG computations 17 , 
was shown to be fairly good in Rcf. [13]. Hence we may determine (3 with a reasonable degree of accuracy from (5.9) 
as long A2 A ' Al <C 1. It then follows that f3 < -j= and hence the SGM (5.3) is in the attractive regime. 
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A. Spectrum of the SGM 

The SGM (5.3) is integrable and many exact results are available. The spectrum of the SGM depends on the value 
of the coupling constant (3. In the so-called repulsive regime, l/y/2 < (3 < 1, there are only two elementary excitations, 
called soliton and antisoliton. These have a massive relativistic dispersion, 



E(P) = ^/M 2 + v 2 P 2 



(5.14) 



where M is the gap. 

In the regime < (3 < l/y/2 relevant to our discussion, soliton and antisoliton attract and can form bound states, 
known as "breathers". There are 



N - 



\-(3 2 
d 2 



(5.15) 



different types of breathers, where [x] in (5.15) denotes the integer part of x. The breather gaps are given by 

M n = 2Msin(n</2) , n = l,...,Af, (5.16) 

where 

P 2 



1 -(3 2 ■ 



(5.17) 



The number of breathers is a function of the applied magnetic field H. There always is at least one breather. A 
second breather appears above a field _ff , which is determined by the requirement 



V3' 



(5.18) 



This requirement is fulfilled for 



H > H k 1.5M . 



(5.19) 



B. Dynamical Structure Factor 

A basis of eigenstates of the SGM is given by scattering states of solitons, antisolitons and breathers. In order to 
distinguish these we introduce labels B\, E>2, ■ ■ ■ Bj^ 7 s,s. As usual for particles with relativistic dispersion, it is useful 
to introduce a rapidity variable 9 to parameterize energy and momentum 

E s {6) = M cosh (9 , P s {6) = {M/v)smh9 , (5.20) 
E 7s {6) = M cosh 9 , P s {9) = (M/v) sinh 9 , (5.21) 
E Bn (9) = M n cosh 9 , P Bn (9) = {M n /v) sinh 9 . (5.22) 

Two-point functions are expressed in terms of a basis of scattering states of solitons, antisolitons and breathers as 
summarized in Appendix B, see Eqn (B8). After carrying out the double Fourier transform we arrive at the following 
representation for the imaginary part of the retarded two-point function of the operator O for u> > 



S°( W ,g) = f;X; f f^CTj l/£.. e „ (ft ■ ■ ■ 8n)\ 2 5(q Yl M * sinh^/g) 6{u> g M ^ cosh ^) ( 5 - 23 ) 

n=l u ^ ' ' j j 

The form factors of the operators exp(±i/3$/2) in the sine-Gordon model were determined in Ref. [28,29]. Using these 
results we can determine the first few terms of the expansion (5.23) for the transverse spin operators. We have 

S**(u,,-+q) = cLf 2 5( s 2 -M 2 )Q(H-H a ) 

a o I! 
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+Rc- 



\F COS (9 )\ 2 



7T 



S"{u,,-+q) 

do 



sVs 2 - AM 2 
C^nf^ 2 - Ml) 
\F sin (9 )\ 2 



+Ro 



sVs 2 - 4M 2 



(5.24) 



Here C is an overall (dimensionful) constant. The terms proportional to F sln and F cos represent the contributions by 
intermediate states involving one soliton and one antisoliton and 



9 = 2arccosh(s/2M). 



(5.25) 



The 5-function contributions are due to the breather bound states. The soliton-antisoliton form factors are given 

by 28,29 

\F^(9)\ 2 = |<O|sin(J$(O))|020i> + - 12 
9(0) 



£cosh(^) 



F^(9)\ 2 = (0|cos(|$(0))|W + 



9(0) 



esinh(^f) 



(5.26) 



where 9 = 9 2 — 9\ and 



g{9) = isinh6</2 

dt sinh 2 (t[l - ie/ir]) sinh(t[£ - 1]) 



x cxp 



(I 



t sinh2i sinh£i coshi 

The absolute values squared of the breather form factors are 28,29 



(5.27) 



fi = 2 sin 
h 



exp 




* £ dt t 



2ir sin t 



2|.g(-^[l-2£])| i 



cot(TrC) cot(^) 2 



(5.28) 



An important result is that the first bound state B\ is visible only in S vv and does not couple to S xx . We note that 
there are additional contributions in the spectral representations (5.24) at higher energies. For example, there is a 
two-breather B\B\ contribution to S xx at energies above 2M\. 

We plot S aa (uj, -^) as functions of uj/M for several values of the applied magnetic field in Figs 3-4. In order to 
give a visual impression of their spectral weights we have broadened the ^-functions corresponding to the breathers 
by convolution with a Gaussian. 

We first discuss the evolution of S lyy (w, ^-) shown in Fig. 3: as H moves away from H c — m the breather B\ 
splits off from the soliton-antisoliton continuum and very quickly takes over most of the spectral weight. Except for 
a narrow window (in magnetic field) above H c the yy-component of the dynamical structure factor is dominated by 
a coherent single-particle peak. 

The evolution of S xx (uj, -^) is very different as is shown in Fig. 4: as H moves away from H c = m the incoherent 
soliton-antisoliton continuum slowly narrows (on the scale of the field-dependent soliton gap) until it eventually begets 
the second, heavy breather B 2 at H = H n s=a 1.5m. Over a large interval of magnetic fields S XK (u>, ^) is dominated 
by the incoherent soliton-antisoliton continuum. 
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FIG. 3: S yy (io, — ) as a function of w/M for several values of the applied field H. 
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FIG. 4: S' xx (w, ^) as a function of lj/M for several values of the applied field H. 

1. Spectral Weights 

In order to compare the spectral weights located in the coherent breather peaks to the spectral weight associated 
with the soliton-antisoliton continua it is useful to define quantities 



M 2 f 25 
— Jo ' 



dx S xx (xM, — ) 

do 



l B 2 



jyy — 



M 2 
~C 



[ dxsyy(xM,-)= /:;, - /?. 

Jo a o 



(5.29) 



For example, CI yy /M 2 is the spectral weight of the yy-component of the dynamical structure factor at the antifer- 
romagnctic wave number integrated over the frequency interval [0,25M]. It has contributions Fg from the coherent 
breather peak and 1JJ from the soliton antisoliton continuum (there are also contributions to due B\Bi two-breather 
states etc, but their contributions are sublcading). It is important to note that the soliton gap M and the overall 
factor C depend on the applied magnetic field. These dependencies drop out once we consider spectral weight ratios 
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such as 

jyy p™ 7"% x 

j-yy ' jyy > ryy ■ l"""^ 

J Bi i Si 1 B 1 

These ratios are plotted as functions of (3 in Fig. 5. We see that for small (that is at H 3> H c ) most of the spectral 
weight is situated in the coherent peak associated with the first breather B\. Very close to the transition the second 
breather does not exist and most of the spectral weight sits in the soliton-antisoliton continua. The crossover between 
these two regimes occurs around (3 ss 0.675, which according to Fig. 2 corresponds to 

— ps 1.025 . (5.31) 

The lesson is that interactions make the summed dynamical structure factor 

S**( u ,l-+ q ) + S yy (w,— +q) (5.32) 
ao ao 



look coherent expect for fields very close to H c . On the other hand, the polarized structure factor S ™^, -£-+q) looks 
incoherent! It would be very interesting to attempt to disentangle the components of the dynamical structure factor 
in inelastic neutron scattering experiments and in this way observe this incoherent scattering continuum. 




FIG. 5: Spectral weight ratios I^f/I^ and Ig*/!^ as functions of the parameter (3. 



2. Polarizations in the LG model 

How do these results fit into the general picture of the LG model? In the latter one expands to quadratic order in 
the fields ipi and ip2 around the minimum of the effective potential at vac = (mo, 0, 0) 

tt2 _ A 2 

ml - ^ . (5.33) 
The effective Lagrangian for the fields (pt t % becomes 

^ 1 / dtp a \ 2 V ( dtfa \ 2 H c dip a 
a—l v / K ' 

H 2 -A 2 2 A2-A? , , s 

— ^ - -hst^ ■ (5 ' 34) 



19 



This is the same as (3.1) for the fields ipi^ if we make the replacement (in (3.1)) 

A 2 — ► 3H 2 - 2A 2 , 

A 2 2 — ► A 2 2 + H 2 - A\ . (5.35) 
This implies for the polarizations in the limit H 3> A2 

A 2 -(0) 



3H Z 



(5.36) 



^i-(O) 
In other words 

|A 2 _|»|Ai_| , (5.37) 

and as a result the coherent low-energy mode is dominantly polarized along the y-direction! This agrees nicely with 
the sine-Gordon calculation, where the dominant feature, the first breather, appears in S yy . 

VI. INTERCHAIN COUPLING 

So far we have considered a purely one-dimensional situation corresponding to an ensemble of uncoupled spin- 
1 chains. As long as the magnon gap is large, a weak coupling between the chains may be neglected in a first 
approximation. On the other hand, the interchain exchange is expected to lead to significant qualitative changes in 
the magnetic response close to the critical point where the magnon gap becomes very small 35 . In order to assess the 
effects of a weak interchain coupling for H w H c wc consider a Landau-Ginzburg model of the form 



C = J2^n+C int , (6.1) 



where 



1 (d(p n - - V v f d<p„ K 

3 A 2 

E o,„2 xi,-? |4 

o=l 

Ant = — y^0j-(fik- (6.2) 

a ° i~Z\ 

Here the sum (jk) is over links between neighbouring sites on different chains and we have dropped quartic terms in 
Lint that arise from the interaction of the smooth components of the spin operators. 

As we have seen in section IV, close to the the Ising critical point the low-energy degrees of freedom are described 
by off-critical Ising models. Hence at low energies we have 

£„ ~ R n d-R n + L n d+L n — irrio R n L n ■ (6-3) 

The leading low-energy projection of the interchain exchange follows from (4.5), (4.6) 



a 



1/4 



where C is a dimcnsionlcss constant. The "quasi-lD Ising model" (6.3), (6.4) has recently been studied in Ref. [36] 
and we may follow some of this analysis here. 

A. Mean-Field Approximation 

As a first step, we analyze the model (6.3), (6.4) by means of a self-consistent mean-field approximation 30,31 . We 
assume the existence of a nonzero expectation value 

0) ^0, (6.5) 
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which corresponds to Neel order along the x-direction. The long-range order can be induced by the magnetic field, the 
interchain coupling or by both. In the presence on a nonzero expectation value (6.5) we may decouple the interaction 
term in (6.4) and arrive at the following mean-field Lagrangian density 

£mf = Rd-R + Ld + L — im RL + — a . 

v 

(6.6) 

Here "magnetic field" h has dimensions of s~ ^ by virtue of the normalization (4.7) and is subject to the self-consistency 
condition 

* (*) , (6-7) 

where Z is the number of neighbouring chains. The mean-field theory is purely one-dimensional and describes an 
off-critical Ising model in an effective magnetic field induced by the neighbouring chains. The model (6.6) has been 
studied by several authors 32-34 and is known to exhibit very interesting physical behaviour as m and h are varied. 
In order to discuss the effects of h and too it is convenient to consider the Euclidean two-point function of Ising order 
parameters 

X* (w, q) = f°° dx dr e^o* (*(t, x)a(0, 0)) . (6.8) 

J — oo 

We note that Xaa ^ related to the xx-component of the staggered susceptibility by analytic continuation to real 
frequencies. The Lagrangian (6.6) defines a one-parameter family of field theories labelled by the dimensionless 
quantity 34 

X = m h~~^ . (6.9) 

In the two special cases \ = and |x| = oo the model(6.6) is integrable and the susceptibility (6.8) can be determined 
to a very high accuracy by means of the formfactor bootstrap approach. In what follows we first review known 
quantitative results for the cases |x| — * oo and x — *• and then summarize the qualitative behaviour for general values 
of X- 



h = ZC 2 v^ 
an 



0() 



1. The Limit h — > 0: McCoy-Wu Scenario 

The regime h — > was studied in Refs 32,33 by means of a pcrturbative expansion in h. In the absence of a field 
(h = 0) the dynamical structure factor has been given in section IV. For mo > the Ising model is in its disordered 
phase and the spin-spin correlation functions are dominated by a single-particle pole and the next-lowest excited 
states occur in the form of a three-particle scattering continuum. The dynamical structure factor is proportional to 
(4.9). Introducing a small magnetic field leads to a small shift in the position of the single-particle pole. Furthermore 
a two-particle scattering continuum of excited states emerges. 

For m < the Ising model is in its ordered phase. This means that there is a nonzero value for the staggered 
magnetization 

(a)l = 2 1 / 6 e- 1 / 4 A 3 TO 74 , (6.10) 

where A denotes Glaisher's constant (4.12). The structure factor in the ordered phase is given by (4.16): the structure 
factor is incoherent and there is a two-particle branch cut starting at uj = 2m . 
It is convenient to define a dimensionless magnetic field h by 

h = { ^h. (6.11) 
ml 

After a resummation of a pcrturbative expansion in h McCoy and Wu established that the spin-spin correlation 
function has the following large-distance behaviour 32 

Wt , iM „, 0)) , w8 ?+™ s 



x ^cxp(-TO r(A i /i) 2 / 3 ) , (6.12) 
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where r 2 = t 2 + x 2 /v 2 and A; are the positive solutions to the equation 

Ji(A ; /3) + J_i(V3) = 0. (6.13) 

The interesting point is that the Fourier transform of (6.12) no longer has a branch cut! There are single-particle 
poles at 

Q 2 +v 2 q 2 = -[2 + {h\ l ) 2 ^] 2 m 2 . (6.14) 

In other words, the two-particle branch cut has disintegrated into a series of single-particle poles. The residues of 
these poles are proportional to 

h [2 + {h\ l f' z ]-\ (6.15) 

Hence the lightest particle carries more spectral weight than the heavier ones. This is quite different from the result 
for h — where the structure factor vanishes as the threshold is approached from above. 



2. The Limit mo — > 0: Magnetic Deformation 

In the limit too — > the model (6.6) is intcgrable 37 . The spectrum consists of eight types of massive self-conjugate 
particles. Three of them have masses below the lowest two-particle threshold. The two-point function of the Ising 
order operator a was calculated by the form factor bootstrap approach in Ref.[38]. The dominant contribution to 
the two point function of Ising order parameter fields is due to the lightest particles. The dynamical susceptibility is 
approximately 



Am 2 



15nh 



2 2 fl 2, (6.16) 



or - Vq* - to' 



where the particle masses are 37,43 



toi w 4.40490858 , 

m 2 w 1.618 mi, m 3 w 1.989 mi. (6.17) 

and 36,38 

Zi w 0.247159 , Z 2 w 0.069017 , Z 3 « 0.02096, (6.18) 

The expectation value of the Ising order parameter is 43 

(a) w 1.07496 h A, (6.19) 

which enables us in principle to solve the self-consistency equation (6.7). The important point is that most of the 

spectral weight is located in the coherent modes corresponding to the two lightest particles. The ratio of weights 
between them is 

^4 - 5.79427. (6.20) 

{Z 2 /m 2 ) 

The region too ~ limit was studied by form factor perturbation theory in Rcf.[34]. 



3. Qualitative Behaviour in the general case 

For general values of \ the qualitative behaviour of Xaai&i l) lii known and may be conveniently summarized 32 ' 34 
by considering the evolution of \aa with \ along a path in the m — h plane as shown in Fig. 6. 

In Fig. 7 we show the analytic structure of Xaa as a function of s = ui 2 + v 2 q 2 for various locations along the path 
set out in Fig. 6. For example, point (a) corresponds to the disordered phase of the off-critical Ising model, where in 
Euclidean space there is a single-particle pole at s = im,Q and a 3-particle branch cut along the positive imaginary 
axis starting at s = 3imo- Point (b) shows the small shift in the position of the pole and the emergence of a 2-particlc 



22 




• * 



2im, 



4im 
2im„ 



FIG. 7: Structure of poles and branch cuts of Xaai^il) at the points (a)-(g) in the h — mo plane indicated in Fig. 6 (mi 
mo[2 + (h\i)i]). 



branchcut 32 . Points (c)-(e) describe the vicinity of the Ising model in a magnetic field; there are several single-particle 
poles below a 2-particle branchcut and the number of these poles increases as we move along the path. Finally, points 
(f)-(g) describe the breakup of the 2-particle branchcut mentioned above. 

An important point is that for mo < 0, which corresponds to the ordered phase of the Ising model for h = 0, the 
general effect of the magnetic field is to make the dynamical susceptibility look more coherent in these sense that the 
low-energy regime is dominated by single-particle poles. In particular, a weak interchain coupling in the ordered phase 
close to H c leads to a disintegration of the 2-particle scattering continuum that dominates the dynamical structure 
factor (4.16) and the formation of a series of single-particle poles. 



B. Beyond Mean-Field: RPA 

It is straightforward to go beyond the mean-field approximation by resumming all diagrams in the interchain 
coupling that do not involve loops. This leads to the RPA expression for the dynamical susceptibility 39 

where Jj_(k) is the Fourier transform of the interchain coupling and 

X xx (w,9) = — [-l 4 X^(w,«). (6-22) 



a 



v 
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In our notations 

Jj_(k) = Jj_[cos(fc y ao) + cos(k z ao)] (6.23) 

for a simple cubic lattice. It was shown in Ref.[36] that in the vicinity of points (c)-(e) of Fig. 6 the RPA leads only 
to slight changes of the mean-field results. More precisely, single-particle excitations corresponding to poles at s = m 
(with residue Z/2m) in the ID susceptibility x xx ( ul ,q) acquire a transverse dispersion ZJ±(k) in the RPA. 

VII. SUMMARY AND DISCUSSION 

We have studied the spectrum and dynamical spin correlations for Haldane-gap systems in the presence of a magnetic 
field. We have paid particular attention to the role played by the crystal field anisotropics present in materials like 
NDMAP. We have concentrated on the case where the magnetic field is applied along the same same direction as the 
largest single- ion anisotropy D (which we identify with the z-direction in spin space). Generalizations of our results 
to other cases is straightforward. Our main results are as follows: 

• At low fields there are three coherent modes M_, M + , M 3 . Their respective gaps A_(_ff) < A + (H) < A 3 (_ff) 
are field-dependent. 

• Above a critical field H d , M 2 develops a finite lifetime via the decay process M + — ► M_M_M_. M_ and M 3 
retain infinite lifetimes. 

• At a critical field H c > the gap A_(_ff) vanishes. In the vicinity of H c the low energy degrees of freedom 
are described by an (off-critical) Ising model and the dynamical structure factor is calculated by exact methods. 
For H > H c the dynamical structure factor is dominated by an incoherent two-particle scattering continuum 
above a finite-energy threshold. 

• At fields sufficiently above H c the low-energy degrees of freedom are described by a sine-Gordon model. 
S yy (uj, ^ + q) is dominated by a coherent single-particle bound state with a spectral gap below a two-particle 
scattering continuum. The most pronounced feature in S xx (lj, ^ + q) is an incoherent two-particle scattering 
continuum above a finite-energy threshold. 

• The effects of interchain coupling are most pronounced in the vicinity of H c . Taking it into account in a mean- 
field fashion leads to a purely one-dimensional effective description at low energies in terms of an Ising model 
in a magnetic Held. This suggests that Haldane-gap materials with single-ion anisotropics in a magnetic field 
may constitute a realization of this very interesting theory. Within the mean-field description the main effect 
of interchain coupling is to generate coherent single-particle modes from the incoherent scattering continua. As 
a result the dynamical structure factor will appear more "coherent" . 

Our findings shed some light on the question why recent inelastic neutron scattering experiments on NDMAP 4,5 
have failed to find any evidence of scattering continua in the high-field phase. At large fields these are suppressed 
through bound-state formation, whereas in the vicinity of the critical field H c the coupling between chains effects 
similar shifts of spectral weight to single-particle modes. It would be interesting to investigate some of our predictions 
experimentally. In particular we hope that it may be possible to 

1. address the issue of the finite lifetime of M 2 above the critical field H&. 

2. disentangle the xx and yy components of the structure factor at high fields. According to our predictions the 
xx-component will remain incoherent up to fairly large fields so that a scattering continuum may be observable. 
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APPENDIX A: LOW-FIELD PHASE IN THE ABSENCE OF CRYSTAL FIELD ANISOTROPY 

In the absence of a crystal field anisotropy the Hamiltonian is 

H(h) = jJ^Sn-Sn+i-hS*. (Al) 

n 

The Heisenberg equations of motion read 

j t St{t) = i[H{%St\TihSt , 

j t S*(t) = i[«(0),S*] . (A2) 

Equations (A2) permit us to express the dynamical susceptibilities for h ^ in terms of the ones in zero field 

X + ~(uJ,q,h) = x + ~(w - h,q,0) , 
X~ + (^,q, h) = x~ + (w + h,q,0) , 

X zz (u>,q,h) = x zz (",q,0)- (A3) 

This implies is that the structure factors in a field are simply the same as in zero field apart from constant shifts in 
energy. The leading contributions to the dynamical susceptibilities in zero field have been calculated in the framework 
of the 0(3) nonlinear-sigma model approximation to the isotropic spin-1 Heisenberg chain in Refs [23,41]. In follows 
from these results that the three-particle contributions are very small. We note that the threshold of the M + M + M_ 
three-particle continuum (two S z — 1 magnons and one S z = — 1 magnon) is at 3 A — h, i.e. for any H < H c it still 
is very slightly higher in energy than the highest energy M_ magnon mode. Hence all three magnon modes remain 
"sharp" for all H < H c . 

We note that analogous considerations apply in the presence of a single-ion anisotropy in z-direction only {E = in 
(1.1) and a magnetic field along the z-direction). The dynamical susceptibilities for finite fields can then be expressed 
in terms of the zero field susceptibilities through the equation of motions for the spin operators. This is quite useful 
for the Majorana fermion model, where the staggered components of the spin operators are expressed in terms of 
Ising order and disorder operators. The latter transform nontrivially under the Bogoliubov transformation used to 
diagonalize the Hamiltonian for nonzero magnetic fields. 

APPENDIX B: SPECTRAL REPRESENTATION OF CORRELATION FUNCTIONS 

In this appendix we collect useful formulas for spectral representations of correlation functions in massive, integrable, 
relativistic quantum field theories. 

We parametrize energy and momentum of single particle states in terms of a rapidity variable 9 

EJ6) = A e cosh 9 , PJ6) = — sinh6» . (Bl) 

v 

Here the index e labels the different types of particles and A e are the corresponding spectral gaps. A scattering state 
of N particles with rapidities {9j} and indices {ej} is denoted by 

\9l, 92, ■ • • , ^A r )ei,e 2 ,---eN ■ 

(B2) 

Its energy and momentum are 

N N A 

mo 3 }) - E A ^ cosh ^ ' p Wi}) = E -f- sinh ^ • ( B3 ) 

j=i i=i 
A basis of states is most easily constructed in terms of the generators of the so-called Faddeev-Zamolodchikov algebra 

^(01)4,(02) - z\ l {e 2 )z\ l {e l )sif 2 {e 1 -e 2 ), 

± z G 2 1 

Z^(9 1 )Zl(9 2 ) = Zl,(9 2 )sf 2 2(9 2 -9 1 )Z^(9 1 ) 

+2tt 525(01-02). (B4) 
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Here S e l' e J(0) is the factorizable two-particle scattering matrix of the integrable quantum field theory. Using the ZF 

€ l' e 2 

operators a Fock space of states can be constructed as follows. The vacuum is defined by 

Z £i (0)\0)=0. (B5) 
Multiparticle states are then obtained by acting with strings of creation operators Z\{6) on the vacuum 

\9 n ...9 1 ) €n ... ei =Zl(6 n )...Z$ i (6 1 )\0). (B6) 
The resolution of the identity in the normalization implied by (B4) is given by 

00 1 rco " in 

1 = EA £/ n?l^---'^.- ei '- en (0i,---,0 n \. (B7) 
n= o n - {£j} °°i=i Zn 

The two point function of some operator O can now be expressed in the spectral representation as 

°° 1 r°° ™ rlf) 

(Ot( t , 3)0(0, 0)) = J2^X II af l/°... eB (fi. ■ • ■ - M 2 exp (-^({^}) + «P({^})), (B8) 

n=0 ' ■'-° j=l 

where the formfactors are given by 

/£..«„ (0i • • • *») = (0|O(0, O)|0„ . . . 0i) e „... ei . (B9) 

APPENDIX C: BOUND STATES IN THE MAJORANA MODEL 

In this appendix we address the question of whether the current-current interaction in the Majorana model leads 
to the formation of bound states. For simplicity we consider only the SU(2) symmetric case with Hamiltonian 

i f 3 

H = - dx^2v[L a d x L a - R a d x R a ] - 2mR a L a 
J 0=1 

+g ( dx^J a J a . (CI) 

J a 

We aim to establish that bound states exist for any g > 0, whereas there are no bound states for g < 0. We recall 
that the Majorana fermion model arises from the spin-Hciscnbcrg Hamiltonian with an additional biquadratic term 

-Hbiquad = J S n ■ S„+l — b (S„ • S„ + i) , (C2) 

n 

where \b — 1| <SC 1. For b > 1 the model is in a dimerized phase whereas b < 1 corresponds to a Haldane spin- liquid 
regime. One may establish by using the expressions (2.5) for the spin operators that the case b < 1 (b > 1) corresponds 
to g < (g > 0). In order to determine whether the current-current interaction leads to the formation of bound 
states, we first consider the limit of a very anisotropic interaction 



i f 3 

Hani = ^ J dx^v[L a d x L a - R a d x R a ] - 2mR a L a 



0=1 

+g J dx J 3 J 3 . (C3) 

This case can be mapped onto a single massive Majorana fermion plus the massive Thirring model by introducing 
complex Fermi fields by 
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The Hamiltonian density is rewritten as H an i = ^Maj + Hmtm , where 



Hmtm 



dx[L 3 d x L 3 - R 3 d x R 3 



2m 



R3L3], 



iv J dx $ R d x $ R - y{d x y L 



+ 



dx 



m[^ R ^ L + h.c] + 2g^[^ L ¥ R ^ R 



(C5) 



In Eqn (C5) we have dropped a term proportional to J dx[i> R ^>n + as it commutes with the Hamiltonian. It 

is well known that in the massive Thirring model there are breather bound states for g > 0, but no bound states exist 
for g < 0, see e.g. Rcfs[44]. 

A different approach is to use large-N methods. If we consider N species of Majorana fermions rather than three, 
we may decouple the interaction through a bosonic Hubbard-Stratonovich field a. For even N and g > the problems 
maps onto the 0(N/2) massive Gross-Neveu model, which is known to have bosonic bound states in the large-N 
limit 45 . 

APPENDIX D: LOW-ENERGY PROJECTIONS OF THE STAGGERED MAGNETIZATIONS 

In this appendix we give arguments in favour of the identification (4.6) at low energies and in the vicinity of the 
Ising critical point at H — H c . We first consider the LG theory and then the Majorana fermion model. 



1. Landau-Ginzburg Model 

It is instructive to examine the evolution of the amplitudes A aa entering the mode expansions (3.7) of the scalar 
fields (fi a as the magnetic field is increased. We recall that the critical field is H c = A\. In the vicinity of H c we 
parametrize 



H = A 1 -S, 5>0. 



(Dl) 



As we are interested only in low energies we may restrict our attention to the "— " modes. From (3.10) we obtain the 
following expansions in S 



A 2 -(0) ' 



2Ai(Af-Af) 
3A? + A% 

(7Aj + Aj)g 
2A!(3A2 + A2) 



Eqns (D2) imply that close to H c we have 



A 2 _(0) cx w_(0)Ai_(0) . 



(D2) 



(D3) 



As we have seen before, close to H c the x-component of the staggered magnetization ipi couples to M_ with a finite 
amplitude Ai_(0) given by (3.11). Furthermore we have the identification (4.6) 



ifi oc a , 

where a is the Ising order parameter field. Eqns (D3) and (D4) together suggest that 

¥2 oc d t cr . 



(D4) 



(D5) 



This claim may be substantiated further in the limit where one of the zero field gaps is much smaller than the other, 
i.e. Ai <C A 2 . As we are interested in energies that are small compared to A 2 , we may "integrate out" the high- 
energy degrees of freedom corresponding to ip 2 in the path integral expression for the staggered magnetization n v . 
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Because H c = Ai is small, we furthermore may take the magnetic field into account perturbatively. The staggered 
magnetization in y direction is 



Averaging n v (t, x) over ip 2 , we obtain 

(n y (t,x)) 2 



n v (t,x) = ip 2 (t,x) . 

^ e iS 2 -2i(H/v) J dt 1 dxi[ip 2 d tl ipi] 



where 



dtdx 



2^ 2 



(dt<P2) ~ 77 (3r¥>2) 



2v v T ' 2 

The leading contribution occurs in first order in the magnetic field 

(n y (t,x)) 2 



~ J dxidti (Ttp2(t,x)(p 2 (ti,x 1 )) 2 dt<pi(h,xi) 

2H f 

= / dx\dt\ G 2 (t - ti,x - xi) d t ip\(t\,x-i) 

v J 



2H 

v 



dx' x d£ x G 2 (i! x ,x\) 



d t (fii(t,x) . 



(D6) 



(D7) 



(D8) 



(D9) 



In the last line we have used that the leading contribution to the integral comes from the region t\ w t, x\ w x. This 
shows that the mixing induced by the magnetic field generates a contribution to n v (t, x) proportional to d t fi(t, x) at 
low energies. 



2. Majorana Fermion Model 

Analogous calculations can be performed in the framework of the Majorana fermion model in the case g a — 0, i.e. 
in the absence of the current-current interactions. In particular, let us consider the case where one of the zero field 
gaps is much smaller than the other 

Ai < A 2 . (D10) 



As the critical field H c = ^/AiA 2 is much smaller than A 2 we may treat the magnetic field term perturbatively. 
We may derive an effective action for R\, L\ only by integrating out R 2 and L 2 (we recall that the third Majorana 
decouples in the absence of interactions) 



Ses ~ Si - -{S 2 H ) 2 , 



S H = M J dx dT[L!L 2 + RiR 2 ] , (Dll) 
where () 2 denotes the expectation value with respect to the second Majorana fermion and 

Si = j ' dr dx[Rid-Ri+Lid+Li-iAiRiLi], 

d ± = "-^f^ . (D12) 
The Matsubara Green's functions are defined as e.g. 

Grr(t, x) = -(T T R( T , x) R(0)) . (D13) 
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Their Fourier transforms are 

iu) + vq 



Gn 2 R 2 (uj,q) = - 
G L2L2 (uJ,q) - 



uj 2 + v 2 q 2 + A 2 ' 
ioj — vq 



2 + v 2 q 2 + A 2 ' 



«A 2 

G R2 L 2 {u,q) = H 2 2 i A 2 • ( D14 ) 

lj z + v z q z + A% 

A straightforward calculation then gives 

C cS = Rxd-Rx + Lid+Li - i[Ai - . (D15) 

In other words, integrating out the second Majorana leads to a rcnormalization of the mass of the first Majorana. 
We note that the dispersion relation that follows from (D15) agrees with the expansion of w_(g) (2.20) in the case 
H <C A2 as it must. What we have shown is that in the case Ai <C A2 it is simply the first Majorana that becomes 
critical at H c . The staggered magnetization in x-direction is expressed at low energies by averaging (2.6) with respect 
to the second and third Majoranas, which gives 

n x {x) cx a 1 {x){n 2 {x)){^{x)). (D16) 

The determination of the operator content of n v (x) at low energies is significantly more involved. In order to obtain 
the low-energy projection, we need to average with respect to the second and third Majoranas 

n»(0,0)«-iir</i 3 > 

x J dr dx{{L 2 (r,x) a 2 (0)) 2 Li(t,x) /xi(0) 

+(R 2 (t,x) a 2 (0)} 2 Rx{t,x) mi(0)} . (D17) 

The expectation values () 2 with respect to the second Majorana can be evaluated by the form factor bootstrap 
approach by utilizing the results of Refs [20-22]. We obtain 



(L 2 (t,x) <7 2 (0)) 2 w De-/4 1 



\/vt + ix 

{R 2 (t,x) <7 2 (0)) 2 w D e-™/*—L=e- mr , (D18) 

where D = TO 1/8 (47r) _1/2 2 1 / 12 e _1 / 8 A 3 / 2 and r 2 = t 2 +x 2 /w 2 . Using these results in (D17) we see that the integral is 
dominated with exponential accuracy by the region r w 0, a: w 0. Hence the operator content of n v is determined by 
the fusion of the disorder operator fi\ with the left and right moving fermions L\,R\. The relevant operator product 
expansions can be worked out following Ref. [46] 

L(t,x)h(0) w -t=^(°) ~ — V^ct(0) + — v^5_ct(0), 



i?(r, x)n(0) w -t=<t(0) - — x/i<r(0) + — V^a+<r(0), 



(D19) 



where z = vr + ix, 7 = exp(— z7r/4)/\/47i : and d T — \{d T =fivd x ). Combining (D18) with (D19) we obtain the desired 
result 

n y cx d T G . (D20) 



APPENDIX E: DERIVATION OF THE SINE-GORDON MODEL IN THE HIGH-FIELD PHASE FOR 

WEAK ANISOTROPY 



In this appendix we show how the Sine-Gordon Hamiltonian emerges as the low-energy effective theory at H > H c 
in the small anisotropy limit A 2 — Ai <C H — H c . We first present a derivation in the framework of the nonlinear 
sigma model and then within the Majorana fermion model. 
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1. Nonlinear sigma model 

The isotropic spin-S Hcisenberg chain in a magnetic field can be mapped onto the 0(3) nonlinear sigma model in 
the continuum limit. Exploiting the integrability of the nonlinear sigma model it was shown in Ref. [13] that for 
H > A the low-energy regime is described in terms of a free boson 

H = ikj dx ^ dx ® )2 + (9x6)2 ] • (E1) 

Here is the field dual to $ and fulfils 

vd x Q = -id T <5> , d T Q = ivd x <f> . (E2) 
The low-energy behaviour of spin correlations follows from the correspondence 

St (-1)" Acxp(±z^e) . (E3) 

The parameters v and (3 were calculated in Ref. [13]. Adding a very small crystal field anisotropy to the Hamiltonian 

Ej2i(S?) 2 - (SI) 2 }, (E4) 

3 

generates a term proportional to 

J rfxcos(/39) . (E5) 

The resulting theory is the sine-Gordon model (5.3). The term D J2j(Sj) 2 merely leads to a small change in (3 which 
we ignore here. 



2. Majorana fermion model 



Our starting point is the Hamiltonian (2.14) describing the two Majorana fermions that couple to the magnetic 
field in the limit A 2 — Ai = to, i.e. vanishing gap anisotropy A = 0. Using 



R,L 



u,n ' ikx (1 \ 

^ e c R . L {k) , 



we may express the Hamiltonian as 



where 



Hi 



A=0 



r (ft,, „ 



M = 



vk + H —im 
im —vk + H 



Now we may carry out a Bogoliubov transformation 



with 



to diagonalize the Hamiltonian. We find 

Wl2 





f dk 


A=0 


1 2ir 



cos(ip k ) -ism(ip k ) \ ( c R (k) 
-isin(^ fe ) cos((p k ) J \ c L (k) 



vn 

tan(2 W ) = — 



(H + sgn(k)\/m 2 + v 2 k 2 ) a\a,k 



+{H - sgn(fc)\/TO 2 + v 2 k 2 ) b\b k 



(E6) 



(E7) 



(E8) 



(E9) 



(E10) 



(Ell) 
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Introducing fcrmions c and d by 



we may express the Hamiltonian (Ell) as 



A=0 



c(k) = a k 8(k) + b k 6(-k) , 
d(k) = b k 9(k) + a k 9(-k) , 



= I 7T \(H + Vm 2 + v 2 k 2 ) c<(k)c(k) 

J 27T L 

+ (H - Vm 2 +v 2 k 2 ) dHk)d(k) 



(E12) 



(E13) 



The low-energy modes occur in the lower band in the vicinity ol ±kp — ±\/ (H 2 — m 2 )/v 2 . They can be combined 
into left and right moving Fermi fields by 



d(x) — cxp(—ikpx)R(x) + exp(ikpx)L(x) . 
The low-energy effective Hamiltonian is then 

H' = iv J dx [tfd x L - F&d x R] , 

where v — v 2 kp/H. We now bosonize the low-energy Hamiltonian using 



(E14) 



(E15) 



L\x) 



1 



2tt 
1 



exp i 



exp 



$(x) + e(x) 

. $(x) - 8(x) 
2\/2 



(E16) 



Here ip and ip arc chiral Bosc fields fulfilling the commutation relations [y(x), f(y)] = 27r«. In terms of the canonical 
Bose field $ = </? + (p and the dual field 9 — ip — (p we find 

H ' = ikj dx [ ( ^ $)2 + {dxQ) ^ ■ (E17) 

The high-energy cutoff in this construction is given by the depth of the Fermi sea in the lower band of (Ell), which 
is H — m = H — H c . So far we have neglected the gap anisotropy, i.e. the term 



^pair — ^A 



^ R ^[ - h.c. 



(E18) 



in the Hamiltonian (2.14). In the next step we take it into account under the assumption that it A is small compared 
to the cutoff H — H c . In this limit 7i pa ir is expressed in terms of the modes as 



. . . dk 



Uk)c\(-k) - c L (-k)c R {k) 



(E19) 



After the Bogoliubov transformation this becomes 

r°° (ik 

-iA / — cos(2<p fc ) [d(k)d(-k] - c(k)c(-k) - h.c] 

JO 27T 



+mixcd terms. 



(E20) 



Dropping the "high-energy" filled band as well as the mixed terms (they contribute in higher orders of A/(H — H c )) 
we have 



Wpair - -*A 



dk 



vk 



27r ^m 2 + v 2 k 2 



[d(k)d(-k) - h.c. 



(E21) 
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Expanding around ifc^ this can be rewritten in terms of the left and right moving fermions as 



H' paXr ^iA^l-^jdx [RL LW] . 



Finally, bosonization gives 



(E22) 



(E23) 



By combining Eqns (E17) and (E23) we see that in the absence of interactions the Majorana fermion model gives rise 
to a sine-Gordon effective theory at low energies in the parameter regime we have been discussing. The parameter 
in (5.3) takes the special free-fermionic value (3 = l/\/2- 

Interactions can be treated in a way analogous to the pairing term. If we drop the interaction terms involving the 
third Majorana (which we assume to have the largest gap) the interaction Hamiltonian reads 



Hint - 2g 3 J 



dx L1L2R1R2 , 



(E24) 



where 53 < 0. Expressing this in terms of the complex fields ^ r^l, carrying out a mode expansion and subsequent 
Bogoliubov transformation and finally projecting to the low-energy band we obtain 



Hint — 2 53 



v 2 kl r 

m 2 + v 2 k 2 F J 



(E25) 



where R and L have been introduced in (E14). Bosonization then gives 

H[ nt ^^Jdx [(d x $f - (d x G) 2 } , (E26) 
where g' 3 = ^j-. This term may be combined with (E17) by rescaling the scalar fields in the standard way 



1 - 93 



$ , e 



1 + .93 



1-53 



e . 



(E27) 



In terms of the rescaled fields the total Hamiltonian takes the form of a sine-Gordon model (5.3) (with a slightly 
changed velocity) where 



1+53 

1-53 



1 1 

V2 V2 



(E28) 



Hence the sine-Gordon model is in the attractive regime. 



APPENDIX F: CORRELATION AMPLITUDE IN THE COMMENSURATE-INCOMMENSURATE 

TRANSITION 



Let us consider the LG model (3.1) in the U(l) symmetric case Ai = A2 = A3 = m. For H > H c = m the low- 
energy degrees of freedom are described by a Luttinger liquid, which can be derived by means of Haldane's harmonic 
fluid approach 25 as follows. Forming a complex Bose field out of the two components of the LG field that couple to 
the magnetic field 



and then bosonizing using' 



25 



*B = 



<£i + i<P2 
V2 



(Fl) 



E 



Am® 



(F2) 
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one obtains, after rescaling the fields <& and 9, a Lagrangian density of the form 



14 



c = 



16tt 



1 fd& 
d \~dt 



dx 



(F3) 



Here v = 2vJ2 



> H—m 



po is the (dimensionless) boson density, which corresponds to the magnetization per site and 

IT is the momentum conjugate to <&. As the LG fields <p a are the staggered components of the spin operators we 
conclude that 



S±oc(-l)^exp(±^) 



(F4) 



where (3 depends on the magnetization and is related to the parameter r\ of Rcf. [14] by 1 — rj. By virtue of (F2) 
the amplitude A is proportional to 

el 

A oc y/p^a ^ , (F5) 

where a is a short-distance cutoff and vertex operators are normalized according to (5.6). The short-distance cutoff is 

a 



a = 



Po 



(F6) 



We note that the short-distance cutoff (F6) diverges as H approaches m as p — > and (F4) describes the asymptotic 
behaviour of spin correlation functions at distances much larger than a. Combining (F6) and (F5) we find that in 
general we have 25 



Act p 2 

Let us now specialize to magnetic fields very close to the critical field H c = in 

H — m -c m. 

As shown in Ref. [25], the parameter (3 tends to so that 

z9 



X( " 1,J ' ' 



where the density is given by 1 



po = —\/2m(H — m). 

7TV 



Combining (F10), (F6) and (F5) we obtain 



A = A'a 4 



H-H c 
J 



(F7) 
(F8) 

(F9) 

(F10) 

(Fll) 



where A' is a numerical, field independent constant and where we have used that oc J. The field dependence of 

the correlation amplitude (Fll) is a universal feature of the C-IC transition. 

Let us apply these ideas to another example of the C-IC transition: the spin-1/2 Heisenberg XXZ chain in a 
longitudinal magnetic field 



n 



(F12) 



where —1<S<1. The model (F12) has a phase transition from a gapless, incommensurate Luttinger liquid phase 
to a gapped, commensurate, spin-polarized phase at a critical value 



H C = J(1 + S). 



(F13) 
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Slightly below this transition, i.e. 

< 1 - « 1 , (F14) 
the transverse correlation functions exhibit the following large-distance asymptotics 

{Sm ^ +1 (0) ~ (-D^ ([gc -J ]/J)i , (F15) 

where A is Glaisher's constant (4.12). This result is obtained as follows: the dependence on the magnetic field is 
universal and given by (Fll). The numerical coefficient is fixed by noting that the numerical results of Ref.[18] show 
that A' is independent of the value of the anisotropy S. Finally, we use that the correlation amplitude has been 
calculated for the free fermion case in Rcf.[47]. The result (F15) is in good agreement with the numerical results of 
Rcf. [18] in the close proximity of the transition, as was already noted in 18 . 
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